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LONG-TERM  GOALS 

In  the  context  of  STRATAFORM  and  EuroSTRATAFORM,  the  long-tenn  goal  of  this  project  is  to 
contribute  toward  simultaneously  understanding  (i)  short  term  oceanic  processes  that  erode,  transport 
and  deposit  sediment  in  the  margin  system  and  (ii)  the  creation  of  the  preserved  stratigraphic 
architecture,  seafloor  morphology  and  sediment  facies  on  continental  margins.  In  order  for  models  at 
such  disparate  time-scales  to  interact,  they  must  communicate  through  expressions  which  upscale  the 
underlying  physical  processes.  Under  Modeling  Task  D5,  the  EuroSTRATAFORM  white  paper 
specifically  states:  “Coherent  techniques  will  be  developed  for  upscaling  individual  processes/events 
into  long-tenn  stratigraphic-architecture  and  seascape-evolution  models.” 

Another  conclusion  common  to  recent  ONR  Coastal  Geosciences  programs  is  recognition  of  the 
dominant  role  played  by  episodic,  high  energy  events  in  driving  both  sediment  transport  and  bed 
formation.  During  STRATAFORM,  the  majority  of  across-shelf  sediment  flux  was  found  to  be 
associated  with  a  few  major  flood/storm  events  occurring  over  just  one  or  two  weeks  every  few  years. 
In  its  section  on  modeling,  the  EuroSTRATAFORM  white  paper  similarly  states,  “Of  special  concern 
will  be  important  events  of  strata  formation  (e.g.,  debris  flows,  extreme  floods)  that  are  difficult  to 
observe.”  Thus  another  long-term  goal  of  this  project  is  to  specifically  understand  the  role  of  energetic 
sediment  transport  in  depositing  sediment  on  margins  and  shaping  morphology. 

OBJECTIVES 

Our  Objective  1  is  to  upscale  results  of  realistic  turbulence  closure  models  to  produce  simple  analytical 
expressions  for  energetic  sediment  transport.  A  wave-resolving  (i.e.,  very  short  time-step),  one¬ 
dimensional  (highly  resolved  in  z)  numerical  model  with  advanced  turbulence  closure  is  being  used  to 
simulate  energetic  suspensions  of  both  mud  and  sand.  The  numerical  simulations  are  being  used  to 
extend  our  recently  derived  analytical  expressions  for  fluid  mud  suspensions  towards  coarser  sediment 
and  reconcile  expressions  for  highly  energetic  fluid  mud  and  plane  bed  sand  suspensions.  For  both  mud 
and  sand  it  is  anticipated  that  strong  damping  of  turbulence  by  stable  stratification  will  play  a  critical 
role  in  determining  the  maximum  possible  sediment  load  for  a  given  energy  level. 

Our  Objective  2  is  to,  in  turn,  upscale  analytical  relations  for  energetic  sediment  transport  to  derive 
simple  expressions  for  equilibrium  shelf  profiles  and  associated  deposits.  Initial  results  of  upscaling 
wave-supported  gravity  flows  indicate  that  the  slope  of  resulting  clinoforms  increase  with  greater 
sediment  supply  and  smaller  wave  height.  For  sandier  environments,  it  is  anticipated  that  onshore 
transport  by  wave  asymmetry  will  balance  offshore  transport  by  gravity.  Preliminary  results  of 
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upscaling  sand  transport  suggest  that  concave  upward  equilibrium  shorefaces  become  steeper  as  wave 
height  decreases,  grain  size  increases  or  wave  period  increases.  The  above  results  of  upscaling  transport 
relations  to  derive  shelf  profiles  are  being  compared  to  global  data  bases  of  shelf  morphology  and  also 
to  the  output  of  other  numerical  models  of  margin  evolution. 


APPROACH 

Objective  1 :  To  provide  guidance  in  upscaling  realistic  turbulence  closure  models  under  energetic 
sediment  transport  conditions,  we  are  running  numerical  model  simulations  of  sediment  stratified 
bottom  boundary  layers  using  the  1-D  General  Ocean  Turbulence  Model  (GOTM,  Burchard,  2002). 
GOTM  is  an  open  source,  FORTRAN-based,  two-equation  (k-epsilon)  turbulence  closure  model 
which  uses  the  most  recent  formulations  for  stability  functions  shown  to  perform  well  in  the  presence 
of  strong  thennohaline  stratification.  We  are  adding  suspended  sediment  to  this  model  and  are 
investigating  model  behavior  as  one  approaches  the  limit  of  critical  sediment- induced  stratification.  We 
are  systematically  examining  the  parameter  space  of  varying  sediment  settling  velocity  (ws)  from 
hindered  settling  in  fluid  mud  (ws  <0.1  mm/s)  to  sand  (ws  >  10  mm/s)  and  the  wave  period  of  the 
bottom  boundary  layer  from  stonn  waves  (order  10  s)  to  tides  (order  10  hours).  We  are  comparing 
these  simulations  to  available  field  and  laboratory  data  of  very  high  concentration  fluid  mud  layers 
within  momentum  deficit  layers  under  waves  (e.g.,  Eel  shelf,  Traykovski  et  ah,  2000)  and  tides  (e.g., 
Amazon  shelf,  Trowbridge  and  Kineke,  1994)  and  of  critically  stratified  constant  stress  layers  of  sand 
under  waves  (e.g.,  Hanover  Wave  Flume,  Dohmen- Janssen  and  Hanes,  2002)  and  of  flocculated  mud 
under  tides  (e.g.,  York  River,  Friedrichs  et  al.,  2000). 

Objective  2:  To  date,  our  upscaling  of  analytical  solutions  to  morphodynamics  relations  for  shelf 
profiles  have  been  based  on  the  following  two  relations  (Wright  et  al.,  2001;  Scully  et  al.,  2002): 

a  g  s  C  p^1  =  cd  lul  ug,  Ri  =  g  s  C  ps"’  lul  2  (1,2) 

In  Eq.  (1)  and  (2),  a  is  the  sine  of  the  bed  slope,  g  is  the  acceleration  of  gravity,  ps  is  the  density  of 
siliceous  sediment  and  s  is  its  submerged  weight  relative  to  sea  water,  C  is  the  depth-integrated  mass 
concentration  of  suspended  sediment  within  the  wave  boundary  layer,  Cd  is  the  bottom  drag  coefficient, 
and  Ri  is  the  gradient  Richardson  number.  The  key  velocities  associated  with  Eqs.  (1,2)  are  the 
downslope  velocity  of  the  gravity  current  (ua),  the  nns  amplitude  of  wave  orbital  velocity  (uw)  and  the 
absolute  amplitude  of  the  instantaneous  velocity,  |u|  ~  (uw“  +  ug  )  ",  all  evaluated  near  the  top  of  the 
wave  boundary  layer. 

For  wave-supported  gravity  flows  on  the  Eel  River  shelf,  Wright  et  al.  (2001)  and  Scully  et  al.  (2002) 
showed  that  observed  transport  and  deposition  rates  could  be  explained  by  Eqs.  (1)  and  (2)  in 
combination  with  a  feedback  mechanism  which  maintains  Ri  in  (2)  near  its  critical  value  of  Ric  =  1/4.  For 
Ri  <  1/4,  turbulence  associated  with  intense  shear  instabilities  suspends  additional  sediment,  increasing  C 
and  Ri,  while  for  Ri  >  1/4,  decreased  generation  of  shear  instabilities  reduces  turbulence  and  causes 
sediment  to  settle.  Because  the  maximum  sediment  load  is  determined  by  the  critical  Richardson  number, 
this  approach  for  predicting  sediment  concentration  is  not  dependent  on  detailed  sediment  or  bed 
properties.  Combining  (1),  (2)  and  the  definition  of  lul  to  eliminate  lul  and  solve  for  ug  and  C  then  yields 
the  following  analytical  solutions  for  the  critically  stratified  case  of  Ri  =  Ric: 

ug  =  uw  (a  Ric  cd_1)  {1  -  (a  Ric  cd-1)2}  1/2 


(3) 


(4) 


C  =  Ric  ps  uw2  g  1  s'1  {1  -  (a  Ric  cd-1)2}  1 

Eqs.  (3)  and  (4)  indicate  that  across-shelf  sediment  flux  due  to  wave-supported  gravity  flows  (ug  C) 
increases  with  wave  orbital  velocity  and  bed  slope. 

WORK  COMPLETED 

In  FY03  we  continued  the  process  of  publishing  our  related  work  from  the  ONR  STRATAFORM 
project  (Wright  et  ah,  2002;  Friedrichs,  2003;  Parsons  et  al,  2003;  Scully  et  al,  2003).  With  respect  to 
Objective  1,  we  successfully  translated  the  GOTM  model  from  FORTRAN  into  MatLab  and 
completed  initial  numerical  experiments  into  energetic,  critically  stratified  bottom  boundary  layers 
(Scully  and  Friedrichs,  2003).  With  respect  to  Objective  2,  we  have  successfully  upscaled  our 
analytical  models  to  derive  relations  for  the  landward  portion  of  equilibrium  clino forms  and  subaqueous 
deltas  (Friedrichs  and  Wright,  2003a, b,c;  Pratson  et  al.,  2003). 

RESULTS 

Our  latest  results  with  respect  to  Objective  2  (Friedrichs  and  Wright,  2003a, b,c)  apply  the  above 
relations  to  solve  for  the  shape  of  the  landward  portion  of  a  stable,  sediment  bypassing  clinoform.  The 
convex  upward  portion  of  a  subaqueous  delta  or  clinoform  subject  to  wave-supported  gravity  flows 
will  be  at  equilibrium  if  there  are  no  across-shelf  gradients  in  gravity  driven  flux  and  the  available  river 
sediment  supply  matches  the  capacity  of  wave-supported  gravity  flows  to  remove  sediment.  In  other 
words,  the  equilibrium  profile  requires 

ug  C  =  Qr ,  (5) 

where  Qr  is  the  supply  of  riverine  sediment  per  unit  distance  along-shelf.  Applying  linear  wave 
theory, 

uw  =  0)  (H/2)  (sinh  kh)'1  ,  (6) 

where  H  is  wave  height,  to  is  radian  frequency,  h  is  water  depth,  and  k  is  wave  number  given  by  the 
dispersion  relation 

k  =  to2  (g  tanh  kh)'1  .  (7) 

Combining  Eqs.  (3)  -  (6)  to  eliminate  ug,  C  and  uw  gives  the  following  relation  for  equilibrium 
bathymetric  slope: 

a  {1  -  (a  Ric  c/)2}'372  =  8  (coH)'3  (sinh  kh)3  Qr  s  g  cd  Ric  2  ps_1 .  (8) 

Eq.  (8)  predicts  that  the  slope  of  an  equilibrium  profile  dominated  by  wave-supported  gravity  flows 
increases  with  greater  water  depth  and  sediment  supply  and  decreases  with  increasing  wave  height  and 
wave  period  (via  k)  (Figure  la).  Equilibrium  slope  increases  with  water  depth  to  compensate  for  the 
effect  of  decreasing  bottom  orbital  velocity.  Slope  increases  with  sediment  supply  simply  because  a 
greater  slope  is  required  transport  a  larger  sediment  supply  offshore.  Slope  decreases  with  increased 
wave  period  because  a  greater  period  decreases  the  decay  of  uwwith  depth.  Finally,  equilibrium  slope 
decreases  with  increased  wave  height  to  compensate  for  the  effect  of  greater  uw.  Figure  lb  displays 
equilibrium  profiles  predicted  by  our  theory  along  with  a  comparison  of  the  theory  to  the  observed  depths 
of  clinoforms  at  STRATAFORM,  EuroSTRATAFORM  and  other  sites.  The  model  offers  a  first-order 
explanation  of  the  delta-front/shelf  profiles  that  exist  off  the  mouths  of  the  Eel,  Waiapu,  Po,  Rhone  and 
Ganges  Rivers  (Friedrichs  and  Wright,  2003a). 


Figure  1.  (a)  Equilibrium  concave-downward  shelf  profiles  predicted  by  Eq.  (8)  (T  is  wave 
period).  The  deepest  profile  is  similar  to  the  Eel,  the  intermediate  profile  is  similar  to  the 
Ganges-Brahmaputra  and  Waiapu,  and  the  shallowest  profile  is  similar  to  the  Po  and  Rhone, 
(b)  1:1  line  along  with  comparison  of  depths  at  which  bathymetric  profiles  reach  a  slope  of 
0.01  (h  01)  as  observed  and  as  predicted  by  Eq.  (8).  From  Friedrichs  and  Wright  (in  press). 


Our  results  with  respect  to  objective  1  are  more  preliminary.  However,  our  initial  implementation  of 
the  GOTM  model  has  been  encouraging.  We  have  successfully  translated  the  GOTM  code  from 
FORTRAN  to  MatLab  and  have  begun  numerical  experiments  including  sediment  induced 
stratification.  Our  initial  results  (Scully  and  Friedrichs,  2003)  suggest  that  several  distinct  types  of 
high  energy,  high  concentration  suspensions  can  be  reproduced  by  the  GOTM  model  with  the  key 
variables  being  sediment  grain  size  and  oscillatory  hydrodynamic  period.  For  cases  where  the  settling 
time  is  long  compared  to  the  wave  period,  fluid  muds  develop  up  to  a  sharp  lutocline  at  the  top  of  the 
oscillatory  boundary  layer.  This  first  scenario  compares  favorably  to  observations  of  hindered  settling 
of  fluid  muds  on  the  Amazon  and  Eel  shelves.  Where  the  settling  time  is  short  compared  to  the 
oscillatory  period,  a  critically  stratified  constant  stress  layer  results,  where  strong  concentration 
gradients  are  present  mainly  near  the  bottom  of  the  boundary  layer.  This  second  scenario  compares 
favorably  to  observations  of  sandy  sheet  flow  under  waves  in  the  nearshore  and  to  suspensions  of 
flocculated  mud  under  tides  in  estuaries. 

IMPACT/APPLICATIONS 

A  present  limitation  in  long-term  modeling  of  continental  margin  evolution  is  realistic  inclusion  of 
hydrodynamic  processes  driving  shelf  deposition.  Based  on  field  observations  collected  over  the  last 
20  years,  complex  wave-averaged  currents  driven  by  winds  and  pressure  gradients  have  been  thought  to 
be  mainly  responsible  for  cross-shelf  sediment  transport  and  flux  convergence  on  energetic  accretionary 
shelves.  Unfortunately,  it  may  be  exceedingly  difficult  to  predict  wind-  and  pressure-driven  near-bed 
currents  with  sufficient  accuracy  to  produce  realistic  deposits  over  geological  time-scales.  The  ONR 
STRATAFORM  project,  however,  recently  identified  a  distinctly  different  mechanism  for  across-shelf 
mud  transport  associated  with  gravity-driven  flows  of  fluid  mud  within  the  wave  boundary  layer. 
Gravity  flows  within  the  WBL  can  be  realistically  modeled  based  on  knowledge  of  fine  sediment 
supply,  approximate  wave  height  and  bathymetry  if  one  assumes  that  the  critical  Richardson  number 
within  the  WBL  determines  the  maximum  capacity  of  the  gravity  flow  to  transport  mud.  Complex, 
externally  forced  mean  currents  do  not  appear  to  play  a  critical  role  in  this  newly  identified  transport 
mechanism.  Thus  the  analytical  model  presented  here  has  the  potential  to  greatly  reduce  the 
complexity  and  computational  limitations  presently  limiting  our  ability  to  perfonn  realistic  long-term 
simulations  of  the  geologic  evolution  of  many  continental  margin  environments. 

TRANSITIONS 

Our  data  on  bed  stresses  and  resulting  sediment  resuspension  from  earlier  years  of  our 
STRATAFORM  project  have  been  made  available  to  others  and  are  being  used  to  verify  various 
bottom  boundary  layer  and  sediment  transport  models.  Our  data  can  easily  be  accessed  via  data 
reports  (which  include  data  summaries  on  diskettes)  and  via  the  VIMS  STRATAFORM  website: 
http://www.vims.edu/physical/projects/CHSD/projects/ONR/index.html.  Published  papers  by  others 
which  have  directly  utilized  VIMS  data  include  Morehead  and  Syvitski  (1999),  Ogston  et  al.  (1999, 
2000),  Reed  et  al.  (1999)  and  Zhang  et  al.  (1999).  Additional  papers  by  non-VIMS  authors 
incorporating  VIMS  data  are  in  preparation.  Our  analytical  formulation  for  sediment  flux  and 
deposition  by  critically-stratified,  gravity  flows  has  already  been  incorporated  into  long-term 
simulations  of  margin  stratigraphic  development  by  James  Syvitski's  group  (Syvitski  et  al.,  2001, 

2002).  Our  analytical  approach  has  also  been  made  available  to  other  modelers,  such  as  Fan,  Harris, 


Niederoda,  Reed,  Swift,  and  Traykovski,  all  of  whom  are  at  various  stages  of  incorporating  gravity 
flows  into  more  complex  numerical  simulations  of  shelf  sedimentation. 

RELATED  PROJECTS 

The  following  active  projects  involving  Friedrichs  also  focus  on  coastal  sediment  transport: 

1.  Forecasting  Scour  Related  Mine  Burial  Using  a  Parameterized  Model.  Office  of  Naval  Research 
(www.vims.edu/physical/projects/CHSD/projects/MBP). 

2.  Sediment  Dynamics  of  a  Microtidal  Partially-Mixed  Estuary.  National  Science  Foundation 
(www.vims.edu/physical/projects/CHSD/projects/CAREER). 
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